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over  the  long-term  case;  however,  the  focal-point  intensity  can  still  be 
several  decibels  down  from  its  vacuum  value,  implying  no  better  than  several 
times  diffraction-limited  performance.  When  the  coherence  length  is  not  much! 
smaller  than  the  diameter,  close  to  diffraction-limited  performance  can  be 
expected . 

Comparisons  have  also  been  made  of  the  reduction  of  on-axis  intensity 
with  no  compensation,  tilt-correction,  and  a full  phase-compensating  aperture 
It  is  shown  that  the  effective  coherence  length  of  the  compensated  aperture 
due  to  the  residual  amplitude  fluctuations  is  greater  than  the  long-term 
coherence  length  by  a factor  of  the  square-root  of  the  Fresnel  number  of  the 
aperture;  for  high  Fresnel  number  systems,  this  larger  coherence  length 
results  in  considerable  increases  in  on-axis  intensity. 


T 


SECURITY  CLASSIFICATION  OF  THIS  P AG  E(H7i  Pn  Data  Entered) 


SUMMARY 


In  this  report,  we  have  calculated  the  degradation  in  the  focal- 
plane  irradiance  distribution  due  to  atmospheric  turbulence,  and  the 
potential  improvement  realizable  by  employing  a wavefront  t ' It-correcting 
aperture.  It  has  been  shown  that,  when  the  aperture  diameter  is  of  the 
order  of  the  outer  scale  of  turbulence,  virtually  no  improvement  is 
realized  relative  to  the  uncompensated  case.  For  the  case  when  the 
long-term  coherence  length  is  small  compared  with  the  diameter,  there 
is  considerable  improvement  over  the  long-term  case;  however,  the  focal- 
point  intensity  can  still  be  several  decibels  down  from  its  vacuum 
value,  implying  no  better  than  several  times  diffraction-limited  per- 
formance. When  the  coherence  length  is  not  much  smaller  than  the 
diameter,  close  to  diffraction-limited  performance  can  be  expected. 

Comparisons  have  also  been  made  of  the  reduction  of  on-axis 
intensity  with  no  compensation,  tilt-correction,  and  a full  phase- 
compensating  aperture.  It  is  shown  that  the  effective  coherence  length 
of  the  compensated  aperture  due  to  the  residual  amplitude  fluctuations 
is  greater  than  the  long-term  coherence  length  by  a factor  of  the 
square-root  of  the  Fresnel  number  of  the  aperture;  for  high  ’’-asnel 
number  systems,  this  larger  coherence  length  results  in  considerable 
increases  in  on-axis  intensity. 
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INTRODUCTION 


To  overcome  the  degradation  In  beam  quality  of  the  far-f ield— or 
focal-plane — intensity  pattern  induced  by  atmospheric  turbulence,  it  is 
necessary  to  measure  and  compensate  for  the  phase  distortions  generated 
in  the  wavefront.  The  simplest  form  of  aperture  compensation  is  the 
correction  for  the  linear — or  tilt — component  of  wavefront  distortion. 
This  component  is  identical  to  the  effective  angle-of-arrival  which 
would  be  induced  in  imaging  the  target  from  the  transmitter  position, 
and  the  concept  has  been  demonstrated,  in  principle,  by  Kerr,  et  at.  [1]. 

The  Improvement  in  beam  quality  so  obtained  represents  the  lowest- 
order  correction  of  a class  of  possible  higher-order  distortions,  up  to 
utilizing  e complete  phase-compensated  or  "adaptive"  aperture  [2,3]. 
Because  of  its  relative  simplicity,  it  is  important  to  quantify  and  com- 
pare the  anticipated  improvement  in  beam  quality  with  that  of  the  full 
adaptive  aperture  for  which,  in  the  presence  of  turbulence  alone  (i.e., 
no  thermal  blooming) , essentially  diffraction-limited  performance  can 
be  expected,  at  least  in  principle.  In  this  report,  we  compute  the 
ensemble-averaged  focal-plane  pattern  for  a gaussian  beam  for  both  an 
uncompensated  and  a tilt-corrected  aperture,  and  the  relative  on-axis 
intensities  for  the  uncompensated,  tilt-corrected,  and  fully  compensated 
apertures. 

In  Sec.  2,  we  calculate  the  tilt-corrected  phase-structure  function 
and  the  MTF  of  the  atmosphere-transmitter  combination.  This  MTF  is 
compared  with  that  derived  by  Fried,  and  it  is  shown  that  Fried’s  MTF 
overestimates  the  atmospheric  transmission  of  high  spatial  frequencies. 
In  Sec.  3,  we  examine  the  effect  of  this  tilt-corrected  MTF  on  the 
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average  ir radiance  pattern  in  the  focal  plane;  in  Sec.  4,  we  derive 
analytic  forms  for  the  M.  for  a gauss ian  beam  propagating  through 
Xolmogorov-type  turbulence.  The  MTFs  are  then  used  to  compute  the 
average  focal-plane  irradiance  as  degraded  by  the  turbulence,  but 
compensated  by  the  tilt-corrected  aperture. 

Of  particular  Interest  is  the  effect  of  retaining  an  outer-scale 
dependence  in  the  turbulence  spectrum.  The  outer  scale  is  the  approxi- 
mate dimension  of  the  refractive-index  correlation  length  in  the  medium. 
Now,  tilt  or  angle-of-arrival  fluctuations  are  produced  by  refraction 
through  the  central  portiop  of  a single  eddy;  hence,  we  expect  negligible 
tilt-component 8 over  dimensions  larger  than  the  outer  scale  of  turbulence. 
The  consequence  of  this  effect  is  that  tilt  compensation  for  aperture 
dimensions  of  the  order  or  greater  than  the  outer  scale  of  turbulence 
produces  insignificant  improvement  in  beam  quality.  This  is  particularly 
important  for  applications  employing  apertures  greater  than  =20  cm  near 
the  ground,  where  the  outer  scale  is  typically  tens  of  centimeters. 

For  the  usually  treated  case  when  the  aperture  is  small  compared 
with  the  outer  scale  cf  turbulence,  we  compare  the  atmospheric  MTF  so 
derived  with  the  approximate  "short-term"  MTF  derived  by  D.  Fried  [4]. 

It  is  shown  that  Fried's  assumptions  result  in  unrealistically  high 
values  of  the  atmospheric  MTF  for  high  spatial  frequencies.  However, 
when  his  approximate  MTF  is  multiplied  by  that  of  the  aperture,  the 
errors  in  Fried's  approach  are  of  the  order  10-percent  optimistic. 

This  results  in  approximately  10-percent  optimistic  results  in  the 
focal-plane  intensity.  Finally,  we  compute  the  focal-plane  patterns 
and  the  ratio  of  the  on-axis  intensity  to  its  vacuum  value  in  terms  of 
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Che  redo  of  che  gaussian  bean  radius  to  the  spherical-wave  coherence 
length  PQ.  The  results  show  that,  even  with  tilt  compensation,  several- 
decibel  loss  can  be  expected  when  wq/po  is  much  greater  than  unity, 
although  the  beam  quality  is  still  better  by  several  decibels  than 
without  compensation. 

To  compare  these  results  with  the  improvement  realizable  from  a 
full  phase-compensating  aperture,  in  Appendix  B we  generate  an  atmo- 
spheric MTF  due  to  the  turbulence-induced  amplitude  fluctuations  which 
remain  after  phase  fluctuations  have  been  compensated.  It  is  shown 
that  amplitude-induced  MTF  has  a coherence  length  larger  than  pq  by  a 
factor  equal  to  the  square  root  of  the  Fresnel  number  of  the  aperture, 
Np.  Hence,  effects  due  to  phase-compensated  apertures  scale  as 
V»oVv  For  the  lar-e  Fresnel  numbers  of  interest  in  high-energy 
laser  applications  for  producing  well  focused  spots,  the  residual 
amplitude  effects  will  be  quite  small.  Specific  examples  are  given  in 
the  text  for  10.6-um  propagation. 
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The  most  general  expression  for  the  phase  corresponding  to  a tilted 
wavefront  Is  given  by  a*jr , where  jr  Is  an  arbitrary  point  in  the 
aperture.  It  can  be  shown  that  2ira/X  - 8 is  the  angle  that  the  plane 
of  the  wavefront  makes  with  the  aperture  plane  (where  X is  the  wavelength) 
The  algorithm  for  choosing  the  vector  a is  that  a*£  represent  the  "best" 
linear  approximation  to  the  phase  over  the  aperture,  and  is  given  in 
Eq.  (9). 

With  ja  appropriately  chosen,  the  tilt-corrected  phase  at  the  aper- 
ture point  _r , $ ' (r) , is  given  by 

$'(r)  ” 4>(r)  - a»£  . (1) 

The  tiit-corrected  phase  difference  between  any  two  points  r^ ,r^  in 
the  aperture  is  then 

<j>'(r2)  - 4 ' (r^)  - <fr(r2)  - 0^)  - a*£  , (2) 

where  £ ■ £2""— 1’  ant*  t^e  tllt_corrected  phase-structure  function  is 
given  by 

d^*(£2,-1)5  0*,^2)_*,^l)l2/> 

* D^Cp)  +<((a*£)2^>-  2<a>*£[iJ>(r2)-<J.(r1)]^> 
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Here, 


y«>  - <t*(i2)-*<£1)]2> 


(4) 


is  the  long-term  phase-structure  function,  and 


Dt(^i*-2) " 2<4*£.^(£2)^^i)]/>"  ^£-)2y 


(5) 


Is  the  change  In  the  phase-structure  function  due  to  tilt  correction. 
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The  effect  of  tilt  correction  on  the  (long-term)  average  Is  to 
"improve"  the  atmospheric  MTF  from  Its  uncompensated  value 

“a*"  " exp[-|D^(r1-r2)]  (6a) 


to  the  tilt-corrected  value 


^ _ 

Fried  [4]  has  added  and  subtracted  2 (a.*£.)2  to  Eq.  (3)  to  write 

it  in  the  form 

D^p)  " ^ (p  ) ~ (£*£.) 2_ 2^.*£. 1 1 ♦ (r 2 ^ “S.*£,2 ^ ^ ^rl ^ 1 ^ 

He  has  then  argued  that  for  any  given  tilt,  as  measured  by  a,  the  devia- 
tions of  the  phase  relative  to  the  tilt,  $'(£)  110  <Kr )“£•£..  is  independent 
of  the  tilt.  From  essentially  this  statement,  the  last  term  in  Eq.  (3*) 
is  set  equal  to  zero.  While  one  can  argue  about  the  reality  of  this 
assumption,  the  most  direct  counterargument  follows  from  actually  cal- 
culating the  quantity  as  given  in  Eq,  (3),  and  comparing  the  result 
with  that  obtained  from  Fried’s  assumption.  A second  powerful  counter- 
argument follows  from  the  demonstration  that  D<J,,  as  calculated  from  the 
first  two  terms  in  Eq.  (3'),  can  actually  be  negative,  which  is  impossi- 
ble from  its  definition.  Fried  has  actually  calculated  a structure 
function  which  can  be  negative,  given  by  the  difference  of  Eqs.  (5.3) 
and  (5.8)  in  Ref.  4j 

D^(p)  - 6.88(p/po)5/3  [l-1.026(p/D)1/3]  (Fried)  , 

but  "approximates"  it  by  his  Eq.  (5.9a)  when  -he  term  1.026  is  set 
equal  to  unity.  We  have  carried  out  precise  numerical  calculations, 
and  shown  that  the  1.026  factor  is  indeed  correct,  indicating  the 
inconsistency  of  the  assumption. 
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For  a spherical  wave  propagating  through  inhomogeneous  (tut  iso- 
tropic) turbulence,  the  phase-correlation  function  is  given  by  the 
expression  [5] 


L • 

B4>)  - - AitVfiy*  ds  J dm  K*n(K,s)  Jo(Kps/L)  , (7a) 

0 0 


and  the  phase-structure  function  by 
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L ® 

D^(p)  - J ds  f dm  K«n(K,8)[l-Jo(Kps/L)]  , (7b) 

0 0 

where  4>n(K,s)  is  the  spectrum  of  refractive-index  fluctuations  at  the 
point  s along  the  path,  L is  the  path  length,  and  k (-2n/X)  is  the  wave 
number  of  the  light.  The  factor  6 varies  from  1 tc  1/2  as  the  range  L 
increases  from  the  near  field  to  the  far  field  of  the  receiver. 

Fried  [A]  has  suggested  that  the  instantaneous  tilt  vector  a_  be 
chose  \ to  jjive  the  best  fit  to  $(r)  in  the  least-square  sense  over  the 
aperture  (at  each  instant  of  time).  For  an  aperture  weighted  by  an 
arbitrary  weighting  function  W(r) , this  implies  that  we  choose  a such 
that 
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For  symmetric  weighting  functions,  it  can  easily  be  shown  that  the 
required  value  of  11  is 


(9) 


In  Appendix  A,  we  substitute  the  expression  for  £ in  Eq.  (9)  into 
Eq.  (5),  and  carry  out  the  indicated  averaging  over  the  turbulence 
fluctuations.  The  result  is 


L » 


Dt  »£.2)  " 8ir2k26  J ds  y*dK  K*n(K,s) 


0 0 


2/W(r)Jl(¥)r2  dr 

_0 . 

00 

J W(r)r2  dr 


tr2-rl 


c°8(r1'r2)1+Jl(-ri)lrrt2 


cos(r2,r1>] 


— 00  — 

CM 

2/w(r)J1(^)r2  dr 

0 

(r2-rl>2 

00 

f n 

J W(r)r3  dr 

r 
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L 

J 

(10) 


Combining  Eq.  (10)  with  the  expression  for  D O^-r^)  given  by  Eq. 

(7b)  and  with  Eq.  (6b)  yields  the  tilt-corrected  mutual  coherence  function 
for  the  field  at  the  aperture  points  ,r^  from  a spherical-wave  source. 
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3.  THE  EFFECT  OF  TILT  CORRECTION  ON  THE  LONG-TERM 
AVERAGE  IRRADIANCE  PATTERN 

We  now  interpret  the  phase  $ of  the  preceding  section  as  <f> (r.£)~ 
the  turbulence-induced  phase  shift  along  the  path  from  the  transmitter 
point  (r,0)  to  the  target  point  at  range  z,  (j>,z).  Then,  neglecting 
amplitude  fluctuations,  the  total  field  at  the  tarp^t  point  is  given 
by  the  extended  Huygens-Fresnel  Principle  [6]  as 

Ufc,.)  * - £ J„a(£)  eik<J>-r)2/2z  d2r  . (11) 


The  quantity  UA(r)  is  the  field  in  the  transmitting  aperture,  and 
k ■ 2ir/X  is  the  wave  number  of  the  radiation.  For  example,  if  we  assume 
the  aperture  field  to  be  a nontruncated  gaussian  of  transmitted  power  P, 
e 1 intensity  radius  wq,  and  focused  at  range  f,  then  the  aperture  field 
is 


(12) 


The  focal-plane  or  far-field  intensity  pattern  is  then  given  by 


I(£,f)  " |u(£,f) 


-±-N2 
B 2 WF 
irw 
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if 


-r2/2w2  --y-£/r  i<p  (j>»r)  2 

e e e dr 


(13) 


where  £ is  the  transverse  distance  of  the  target  point  from  the  optic 


axis,  and  N_  ■ kw  /£  is  the  Fresnel  number  of  the  aperture. 
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In  the  absence  of  turbulence,  there  are  no  phase  fluctuations — 

$ ■ 0 — and  Eq.  (13)  produces  the  vacuum  intensity  in  the  focal  plane: 


P 2 
2 WF 

TTW 

O 


-p2/(wo/Np)2 


(14) 


Let  us  now  assume  that  the  dominant  phase  distortion  induced  over  the 
plane  of  the  transmitter  t_,  from  a point  source  located  at  the  focal 
point,  is  a pure  tilt  of  the  wavefront,  or  angle-of-arrival  fluctua- 
tions. It  can  then  be  seen  that 
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♦t  ■ k(0‘r) 


(15) 
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where  6 is  the  instantaneous  arrival  angle.  If  the  arrival  angle  is 


independent  of  the  location  of  the  point  jg_,  then  substituting  Eq.  (15) 
into  Eq.  (13)  yields 


-(£-«() 2/(wo/Sf)' 

I<£,f)  '-Li'.  0 = I <£-«£. f)  • (16) 
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Comparing  Eq.  (16)  with  Eq.  (13),  it  follows  that  if  1)  amplitude 
fluctuations  are  negligible,  2)  the  phase  distortion  from  a point  in 
the  target  plane  is  pure  wavefront  tilting,  and  3)  the  same  tilt  angle 
would  be  generated  by  all  points  in  the  target  spot  dimension  (i.e., 
the  spot  lies  in  an  isoplanatic  region) , then  the  degradation  takes  *he 
form  of  a diffraction-limited  spot  that  dances  with  time  in  the  focal 
plane. 
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Now,  is  the  above  three  conditions  are  met,  we  can  introduce  a com- 
pensating tilt  in  the  transmitted  wavefront  by  changing  the  aperture 
field  from  U.(r)  to  U.  (r)  e”ik(§‘£)4  The  cancellation  of  the  phase 
terms  in  Eq.  (1)  then  implies  complete  compensation  for  beam  dancing  in 
the  target  plane.  To  compute  the  effect  of  tilt  compensation  for  arbi- 
trary <f>(£,r)  on  the  focal  plane,  or  far-field  pattern,  we  write  Eq.  (13) 
as  a double  integral: 


*<£•»  ■ {stffj f Wte 


) e 


ik  , . 

~ T — 1~— 2 


• e 


i[4>  (£,  r ^ £2>1  2 2 


d _r x d‘r2 


(17) 


The  average  turbulence-degraded  intensity  pattern  is  found  by  averaging 
Eq.  (17)  and  using  the  relation 


\e  / • ■ “4-  2 V-r£2>J  ■ 


(18) 


where  Mg  is  the  spherical-wave  mutual  coherence  function,  and  is  the 
phase-structure  function.  For  the  long-term  average,  is  a function 
of  the  difference  coordinate  £ ■ 2’  substituting  in  Eq.  (17)  and 

letting  R ■ 2^  y*e^8  ^or  the  average  intensity 


ik 

- 1"  £“£, 
e M^) 


(19) 


where  the  total  long-term  MTF  is 


Bk: 


W^.!  nw  'A  ' JUJJHS.JH  UW4JWVU  JU  IWJ^WIISipiRWWfpPPIIi 
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M 


«T  (£> 


- f yp> 


/ua(i+U)u a(r-U) 


ik  „ 
e d R 


(20a) 


Aside  from  a constant  factor,  the  integral  in  Eq.  (17)  may  be  recognized 

.it 


as  the  MTF  of  the  aperture  distribution;  hence,  MT  is  the  product  of 
the  atmospheric  MTF  and  the  MTF  of  the  transmitter.  For  the  tilt- 
compensated  aperture,  the  effective  phase  distortion  at  the  aperture  id 
reduced  from  4>(r)  to  $'(r);  as  shown  in  the  next  section,  the  resulting 
phase- structure  function  is  a function  of  r^  and  r^*  or  both  £ and  R. 
Then,  for  the  tilt-corrected  aperture,  the.  average  intensity  is  again 
given  by  Eq.  (19),  but  with  M^  given  by 


X'  '/ ua(Hf  £■)“*(-- H 


-!v<a.P  rt2„ 

e e r d R 


-1VC> 


/ »A(s+^  1 p) 


Ti'i  f Dt<£>£)  , 

e e a k 


(20b)  | 

S 

A 

In  contrast  with  the  long-term  MTF  case,  because  the  tilt-corrected  > 

•1 

phase-structure  function  is  not  a function  of  the  difference  coordinate  1 

£ alone,  the  total  tilt-corrected  MTF  cannot  be  written  as  the  product  1 

.i 

•1 

of  an  atmospheric  MTF  and  an  aperture  MTF.  However,  the  total  distortion  \ 

l 


can  still  be  calculated  from  Eq.  (20b . 
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4.  THE  GAUSSIAN  BEAM  AND  KOLMOGOROV  TURBULENCE 

4.1  The  Tilt-Corrected  Structure  Function 

For  the  gaussian  beam  of  Eq.  (.2),  the  weighting  function  over  the 
aperture  (aside  from  a normalization  factor)  is 

-r^/w^ 

W(r)  |UA(£)|2  - e 0 (21) 

go 

Substituting  in  Eq.  (10)  and  letting  x ■ ^ yield3  for  the  tilt- 

correction  to  the  structure  function 

Dt (£»£.)  " 87T2k26  j ds  /dx  x.aj  H(x,£,R)  , (22) 

0 ° 0 


where 


H(x,P,R)  - 2f(x)  G(x,£,R)  - f2(x)  (p/wo)2 


(23) 


.Here, 


w 

J W(y)  J1(2xy)y2  dy 


f(x)  - 


W(y)y  dy 


x e 


-x 


for  W(y)  ■ e y , and  (noting  that  r^  - R + £_/?,  m h.  “ Sj 2) 


(24) 


where  is  the  Bessel  function  of  order  1.  We  have  shown  numerically 

that  the  Bessel  function  may  be  approximated  by  the  first  two  terms  in 

* 

the  series,  yielding 


2 2 

H a (p/wQ)2[2xf (x)-f2(x)]  - ~^2+'\  (l+2cos2t|i)  x3f(x)  , C 

4w  v 

LOO  J 


where  ^ is  the  angle  between  £ and  J*,  and  f (x)  is  given  by  Eq.  (24) , 
For  a Kolmogorov  turbulence  spectrum,  we  let 


*n<K,o) 


0.033C 
n 

, ,111/6 

K +K 

o 


where  C is  the  index  structure  constant,  and  K *»  L is  the  spatial 
n o o r 

frequency  corresponding  tc  the  outer  scale  Lq.  Allowing  both  Cn  and  L 
to  vary  along  the  propagation  path  and  substituting  Eqs.  (27)  and  (2^ ; 
into  Eq.  (22)  yields 


Dt(p,r)  - 0.82k26wo5/3(p/w q)2  J C2(s) (s/L)5/3  B(s)  , ( 

0 


where 


• * 2 2 
Under  Fried's  assumption,  H would  be  given  by  ■ f (x)(p/wQ)  . 
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B(Q)  - A1(s)  - A2(b) 


d2  R2  2 

— (1+2cob S) 
4w0  wo 


(29) 


Here, 


A^(tf) 


a 

"i  / 


e~y  (2-e-yl 


o ty^0(s)] 


11/6 


dy  - e 


-(i  ■ yo)  -y0r(-  i > yo) 


2-7/6  e2»o 


r(i>2yo)‘2yor('l’2!,o) 


(30a) 


and 


A0(s) 


a 

u 


?.  -y 

jl‘'s”u76  dy 


0 [yty0(s)3 


2 * 


r(?'yo)-2ycr(iyo)+yor('l-yo) 


(30b) 


where 


w B 
o 


yo  ■y0(E)  ‘ \ 2L_L 


(31) 


and 


T(x,y)  * J e”fc  ty_1  dt 
x 


is  the  incomplete  gamma  function 
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From  Eqs.  (28)  through  (31),  we  see  that  takes  the  form 


(32) 


f 

2 2 ' 

j B1  - B2 

+ — (l+2cosZi|») 
4w  w 

l 

Loo  J 

4 

where 


Li 

B.  , - 0.4lk26w3/3  ( C2(s)(s/L)5/3  A.  ,(s)  ds  . (33) 

1,2  o J n 1,2 

0 

4.2  The  Tilt-Corrected  MTF 

We  now  compute  the  overlap  integral  of  Bq.  (20),  using  Eq.  (32) 
for  Dt  and  assuming  a gausslan  beam.  Thus, 


m£C  - e 


2 


2 

1 eZ  C dZR 


-^D  (p)/«-p2/4w2  -r2/w2 

e I e e exp< 


w 


b2p2 

2 

_£_ 

CNJ 

3 

, 2 
4w 

° 1 

o 

(1+2COS  \J)) 


)d2r 


"ID^(p)  2 ~p2/wo  Blp2/2wo  “B2pA/4wo 

a tiw  e e e 

o 


oo 

f 


-x(l+2B2p2/w2) 
e x dx 


-|d  (p)  -p2/v2  B P2/w2  -B2p4/4w‘ 

e ttw  e e e 

o 
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-7Di(P)  “P 


Note  that  the  term  e 


' 2 
ttw  e 
o 


Is  the  result  that  would  have 


been  obtained  for  the  overlap  integral  MT  (£,)  of  Eq.  (20a).  Hence,  we 
obtain  for  the  effective  tilt-corrected  atmospheric  MTF: 


-|yp>  v2/wo 

. “ * a a 


, ,2  / J 

22f 

•\v  2 w0/ 


where  2 is  given  in  terns  of  A^2  by  Ecl*  (3*) . 


A . 3 Approximate  Solutions 


w « L 
o o 


We  first  examine  the  frequently  considered  case  where  the  aperture 
dimension  is  small  compared  with  the  outer  scale  of  turbulence.  In 
this  case,  the  long-term  structure  function  is  given  by 

R 

D+(p)  - 2.91k26p5/3  J C2(s)(s/L)5/3  ds  £ 2(p/pq)5/3  , (36 

0 

where  P is  the  e_1  point  of  the  long-term  MTF.  Then,  for  w0/LQ  ■+  0, 
o 

Eq.  (30)  yields  “ T(l/6)(l-2  ) " 3.08  and  A2  - “ 0.A6; 

from  Eq.  (33),  B1>2  -*»  0. 28(wq/po)5/3  A^.and  B1>2  - (0.86,0.13)  (wq/po)‘ 
Substituting  in  Eq.  (35)  then  yields 


10/3  —1/ 2 

<°  ■ [i+o-52© ' (t)2+n-o56)  fe). 

. exp  -(p/p0)5/3  l-0.86(p/vo)1/3  + 0.13(p/wo)7/3j  , 

k 
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f 

i 

i 

Ir 


I 


and  the  total  tilt-corrected  MTF  of  the  atmosphere  and  lens  combination 
is  given  by 


2/  2 
-p  /wo 


(38) 


For  comparison,  the  long-term  MTF  is  given  by 


-p2/v2  -(p/p0) 

e e 


5/3 


(39) 


From  Eq.  (39) , we  obtain  the  usual  result  that  turbulence  only  degrades 
the  beam  quality  when  pQ  < wq;  i.e.,  when  the  coherence  length  is  less 
than  the  aperture  dimension.  In  Fig.  1,  we  plot  the  tilt-corrected 
MTF  given  by  Eq.  (37)  for  various  values  of  the  parameter  (p0/wq). 

To  compare  these  results  with  those  that  would  be  obtained  using 
Fried's  assumption,  we  use  the  fact  that  Fried  would  have  obtained 
Eq.  (22)  with  H(x,P^,R)  replaced  by  Hp  - f2(s) (p/wo)2.  Then,  Eqs.  (28) 
and  (29)  remain  valid  with  * C,  and 


w 

kf 


y c"2y  , .1 

v ^ 771 


ra.6)  - 2.48 


Tnen,  ■ 0.28(wo/po)5^  « 0.70(wo/po)^^  and,  hence,  Fried's 


atmospheric  MTF  becomes 


■i 


! 


3 


1 

• 3 


(40) 
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and  Fried's  total  MTF  becomes 


(41) 


In  general,  while  the  tilt-corrected  atmospheric  MTFs  predicted 
by  our  more  precise  theory  are  very  different  from  those  generated  with 
the  Fried  approximation  (as  indicated  in  Fig.  1),  when  we  form  the 
composite  MTF  of  the  aperture  atmosphere  system,  which  is  the  only 
measurable  quantity,  the  results  are  not  greatly  divergent.  Generally, 
Fried's  assumptions  produce  considerable  errors  for  the  higher  spatial 
frequencies  of  the  atmospheric  MTF,  in  the  region  where  the  MTF  of  the 
aperture  1j  small,  which  reduces  the  effective  error.  As  a result, 
the  total  tilt-corrected  MTFs — plotted  in  Fig.  2 — differ  primarily  in 
the  region  for  which  the  transverse  dimension  p 1b  greater  than  the 
gau8slan  radius  of  the  aperture,  wQ. 

We  compute  the  focal-plane  irradiance  patterns  for  the  respective 
cases  by  taking  the  Fourier  transform  of  the  total  MTF.  From  Eq.  (19), 
for  Mt  independent  of  the  direction  of  £,  we  obtain 


<I(a,f)> 

I (oTf) 

vac  ’ 


00 

2 /mt(x^wo/po)  JQ(2ax)x  dx 
0 


(42) 


for  the  average  intensity  normalized  to  the  on-axis  vacuum  pattern.  The 

parameter  a ■ kpw  /2f  ■ irpv  / Xf  is  the  ratio  of  the  angle  subtended  by 
o o 

the  point  p at  the  aperture  to  the  vacuum  diffraction,  angle. 
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In  Fig.  3,  we  plot  the  irradiance  pattern  of  Eq.  (42)  us  in,’  the 
total  tilt-corrected  MTF  and  comparing  the  result  with  Fried’s.  The 
numerical  results  yield  an  approximately  10-percent  difference  in  the 
on-axis  intensities,  with  Fried’s  being  optimistic  due  to  the  higher 
spatial  frequency  content  of  his  MTF.  In  Fig.  4,  we  plot  the  reduction 
in  focal-plane  intensity  relative  to  the  vacuum  value  versus  wQ/p0 — 
the  ratio  of  the  beam  radius  to  the  coherence  length  pQ.  From  Eq.  (42), 
this  reduction  io  given  by 

<1(0, f))  r 

- 2 J Hr(x.wo/PJ)x  dx  • <«! 

vac  * q 

To  estimate  the  improvement  anticipated  with  a full  compensating 
aperture,  we  derive,  in  appendix  B,  an  approximate  MTF  which  contains 
the  effects  of  the  amplitude- induced  distortion  remaining  after  phase 
compensation  has  occurred.  It  is  shown  that,  in  contrast  with  the  phase 
result  where  turbulence  degrades  the  pattern  for  pQ  5 wo,  amplitude 
effects  appear  when  pQ  » wQ/2^Np,  where  Np  is  the  Fresnel  number  of 
the  lens.  For  the  large  Fresnel  numbers  of  interest  in  high-power 
applications,  the  amplitude  effects  are  indeed  much  smaller  than  the 
phase  effects,  indicating  considerable  increase  in  on-axis  intensity 
with  a full  compensating  aperture.  For  illustrative  purposes,  we  have 
chosen  a sma31  Fresnel  number  of  4.75,  corresponding  to  wQ  ■ 20  cm, 

A ■ 30.6  pm,  and  f • 5 km.  The  results  for  all  of  the  cases  are  shown 
in  Fig.  4,  indicating  the  increase  in  on-axis  intensity  for  the  various 
configurations.  It  is  seen  that  the  Fried  solution  considerably  over- 
estimates the  Improvement. 


Relative  Intensity  patter 
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5'.  ? r 


Normalized  off-axis  angle,  a <= 


Fig.  3b--Relative  intensity  patterns  for 
tilt-corrected  apertures  (Fried) 


Fig.  4 — On-axis  intensities  relative  to  vacuum  for  uncompensated, 
tilt- corrected,  and  full -compensated  apertures 
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For  the  case  where  the  outer  scale  of  turbulence  is  large  compared 
with  the  aperture  dimension,  the  long-term  phase-structure  function  can 
be  well  approximated  by 


■ 2<*2> 


r 2,  ,2 

-p  /10 

1 - e 


(44) 


where,  for  uniform  turbulence, 


^2^>«  0.39k2C2L3^3  AL 


is  the  phase  variance  at  a point.  We  then  evaluate  the  functions  A 
of  Eq.  (30)  using  the  asymptotic  expansion 


1.2 


_y 

r(a,yo)  e y6”1  e ° , yo  » 1 

From  Eq.  (30),  we  obtain  ^ “ (7/8)y~11/6,  A£  * y~11/6. 

The  condition  yQ  » 1 follows  from  the  assumption  wq  » Lq, 
assuming  that  the  turbulence  is  contained  not  closer  than  a distance  L 
from  the  focal  plane,  where  L1  i LLq/vq.  Then,  from  £q.  (33), 


11/3  L 


1,2 


(0.36,0.41)k26w3/3^j  j C2(s)(L/s)2  ds 


(45) 


and,  for  Cr  constant  between  and  L, 


..  --'i  iv  ii.  .'2  V'  . 


to  dominate  the  linear-phase  fluctuations  for  apertures  of  these  j 

dimensions,  implying  little  improvement  when  they  are  removed.  j 
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Appendix  A 

CALCULATION  OF  THE  TILT-CORRECTED  STRUCTURE  FUNCTION 

In  this  appendix,  we  evaluate  the  second  and  third  terns  of  Eq.  (3) 
(p.  4)  for  an  arbitrary  symaetric  aperture  weighting  function  W(r). 

Evaluation  of  <^(a,£)2^> 

The  vector  a defined  by  Eq.  (9)  of  the  main  text  lies  in  the 
aperture  plane,  and  for  isotropic  turbulence  will  be  uniformly  distri-  ■ 
buted  in  angle  over  the  plane  of  the  aperture.  For  isotropic  turbulence, 
the  expected  value  ^coG2(a,£_)^  *1/2,  where  (a,£)  denotes  the  angle 

it 

between  £ and  and  we  obtain 

<(<£*£.) 2)>  " p2  . (A-l) 

We  write  Eq.  (9)  as  a double  integral,  which  yields 


m “ tffa  d2I2  W<£.i>V,(i2)  <(+(£1)  <Kr>2^  • (A~2) 


where 


A 


it  J" W(r)r^  dr 
0 


The  term  in  Ecl*  (A“2)  is  the  correlation  function, 

B^^-r^),  between  the  phases  at  the  points  r^  and  r^.  Substituting 

it 

The  same  result  can  be  derived  directly  from  the  definition  of 
£ given  by  Eq.  (9). 
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the  expression  for  a spherical  wave  source  propagating  through 
inhomogeneous  turbulence-r-Eq . (7a)— -and  inverting  the  order  of  inte- 
gration, we  obtain 


L °° 


<•*>-  ds  J dK  K*n(K,s)  /jfd2*!  .2  W(I.i)WC£2)  ^-1-2^ 


0 0 


• jo(t  Is^l)  • 


(A-3) 


into  which  we  substitute  the  identity 


Jo(t-  cos[»(r1(r2)]  J r,)  , (A-4) 

where  eo-l,  en«2,  n*>l,2,...  . The  term  contains  the  term 

cosCr^.r^)  and,  because  of  the  angular  orthogonality,  the  only  nonvanish- 
ing contribution  comes  from  n«l.  Thus,  we  obtain 


L « 


<(a2]>»  4Tr2k26E;LA  J ds  J dK  K$n(K, sjjfJ'd2^  d2r2  WCr^WCr^r^ 


0 0 


cos 


(jLl,X2)  Ji(l  rl^  Jl(lf  *2) 


(A-5) 


Expanding  cos‘J(r^,r2)  in  terms  of  the  double-angle  formula,  only 
the  "1/2"  remains  after  integration,  and  we  obtain 


w w 

■|<a2)>p2  - 2w2k26p2yds  J dK  K*n(K,s) 


0 0 


• 2ttA  / W(r)  J.  rlr  dr 


(A-6) 


for  symmetric  W(r). 


Evaluation  of  ^ ■ 2 (r2)-<f>  (r^) 


With  £ defined  by 


-a/. 


W(r)<J>(r),r  d £ , 


(A-7) 


we  have 


2A ^ d2jr  W(r)  ^(£*  r.2>"(— ‘-r-l)]  <(4> (r.) <t>  (±2 ' ” $ CtH  (Ei 


2A J* i2r  W(r)  [(r  • r^  (r-r2)  + <r -r^B^-r^ 
" (£\Li> \ <iT£2)  ‘ B4  (^1}] 


(A-8) 


Substituting  Eq.  (7a)  for  B ^ into  Eq.  (A-8)  and  interchanging  the  order 
of  integration,  we  obtain 


L » 


- sA25 /ds/dK  K*n(K,s)[r(t1,i1)+r(£2,j:2)-r(£1,i2)-r(£2,x1)  , 


0 0 


(A-9) 
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i 


I 


( 


* 


t 


* 


» 


where 


» 2w 

" Ari/dr  r2  w<r>  y d«  cos^r^)  JQ^  li'^jl)  • 


If  we  substitute  Eq.  (A-4)  in  Eq.  (A-10)  and  expand  the  cosine 
as  cosl^r.rj)-^,^)],  we  obtain 


QO 

T(l±»I j)  “ ZttA^  cos(r1>r^)  j W(r)  dr 


Collecting  terms,  we  thus  have 


ip  - 8Tr2k26 /ds  /dK  K4n(K,s)  • 2itA  ^WCr)  ^(^jr2  dr 
0 0 0 


c08(Ll»£2)]+J 


[rl“r2 


g°  s ( , r_2  > J | 


!>  ft  '( • t-K'X.  tl  J v ir-hiteA 


sX  fcW,! v.' RV.“. 


(A-10) 


term 


(A-ll) 


(A-12) 
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Appendix  B 

THE  PHASE- COMPENSATED  MTF 


With  the  compensation  of  phase  fluctuations,  the  MTF  is  given  by 
the  amplitude  term 

- ? “ --KtxCrJ-xfc)]2)  -<X2)  d-B  (p)] 
M-eZX-e/N  1 1 - e w x 


(B-l) 


where  (r^)^>  ■ (*.2^  “ 1b  the  variance  of  log  amplitude,  and 

is  the  amplitude  conduction  function.  How,  it  is  known  that  for 

<*2>«- 


<X2>  ■ y f)  - 0.13k7/6C2z11/6  • 0.25^z/kp2J/6 

/ 2\S/6  5/3  (wo/po)5/3  ✓ -s 

' 0.25(z/kw2)  (w0/P0)5/3  ■ -^576—  • <X'>  « 1 • 0M> 


where 


(B-3) 


2 

is  the  spherical-wave  coherence  length,  and  N_  ■ kw  /z  is  the  Fresnel 

r O 

number  of  the  lens.  For  ~ 1*  the  variance  of  log  amplitude 

saturates  to  a value  of  order  unity;  hence,  we  have,  as  an  iroximate 
formula. 
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(B-4) 


/ 2\ 

From  Eq.  (B-4),  it  follows  that  <^x  / becomes  appreciable  when 


”0  4 


In  contrast  with  the  phase  effects,  which  become  important  when  pQ  s wQ, 

much  stronger  turbulence  is  requried  before  amplitude  effects  appear. 

The  MTF  will  only  differ  appreciably  from  unity  in  the  saturation 

regime,  when  ^(2^  * 1.  In  this  regime,  it  has  been  shown  [7]  that 

the  correlation  length  of  the  amplitude  fluctuations  is  pQ  itself, 

2 2 2 

and  we  approximate  the  correlation  function  by  PQ/(P  +pQ).  Hence,  an 
approximate  MTF  which  retains  the  essential  dependence  on  the  param- 
eters of  the  problem  is 


M - exp-j 


4n: 


576' 


1 + — 


<wo/po> 


2 . 2 
P + P„ 


573 


- exp<- 


4N 


57T 


1 + 


(wo/po) 


573 


[(wo/po)x]2 
1 + [(wq/p0)x]2 


(B-5) 


where  x ■ p/w  . 

o 


■ ( iv 


'"S',  ■ ■ 

■ •■«  • - - 
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